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Abstract 



We study the first-order finite-temperature electroweak phase transition 
of the SU(2) gauge-Higgs model defined on a 4-dimensional isotropic lattice 
with temporal extension Nt = 2. Finite-size scaling study of Lee- Yang zeros 
yields the value of the Higgs self coupling of the endpoint at A c = 0.00116(16). 
An independent analysis of Binder cumulant gives a consistent value for the 
endpoint. Combined with our zero-temperature measurement of Higgs and 
W boson masses, this leads to Mh, c = 73.3 ± 6.4 GeV for the critical Higgs 
boson mass beyond which the electroweak transition turns into a crossover. 

PACS number(s): ll.10.Wx, 11.15.Ha, 12.15.-y 
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I. INTRODUCTION 



The Minimal Standard Model predicts that the electroweak interaction undergoes a first- 
order phase transition at a finite temperature for light Higgs boson masses. A focus of recent 
studies has been whether the first-order phase transition survives with sufficient strength for 
a realistically heavy Higgs boson mass |]], since the feasibility of electroweak baryogenesis 
depends crucially on it. 

The first-order nature of the electroweak transition for light Higgs bosons can be shown 
within perturbation theory. However, perturbation theory breaks down for Higgs boson 
masses larger than about M w due to bad infrared behavior of the gauge-Higgs part of the 
electroweak theory ||. Hence numerical simulation techniques are needed to analyze the 
nature of the transition for heavy Higgs bosons. 

Extensive studies in this direction have already been performed within the effective 3- 
dimensional theory approach, in which all non-static modes of the system are integrated out 
perturbatively. This approach has the advantage that the full Standard Model including 
fermions can be mapped onto a 3-dimensional SU(2) (or SU(2) (g) U(l)) gauge-Higgs model, 
as there are no fermionic static modes at finite temperature. In addition, thinning out the 
degree of freedom to those of a 3-dimensional theory significantly reduces the computational 
requirement. 

Results from simulations in this approach show that the first-order electroweak transition 
weakens as the Higgs boson mass increases , and that it turns into a continuous crossover 
for heavy Higgs bosons with a mass Mh Mw 0- Detailed studies of the endpoint of the 
first-order transition including its universality class have also been made |]8HlO|. 

A potential problem with the 3-dimensional approach is that it relies on perturbation 
theory to derive the 3-dimensional action so that numerical predictions may involve sys- 
tematic errors due to truncation of perturbative series. From this point of view a direct 
simulation of the 4-dimensional system is preferred. Results from 4-dimensional simulations 
provide a check on those of the 3-dimensional method. 

Early studies of the 4-dimensional SU(2) gauge-Higgs system were carried out in 
Refs. [|TT|-|T5| . More recently advances have been made with the use of space-time anisotropic 
lattice [0,[nj • This approach alleviates the double-scale problem that there are light modes 
with long wave length, £ >> l/T, near the endpoint where the transition is of second order. 

In this article we report on a study of the endpoint of the SU(2) gauge-Higgs model 
employing 4-dimensional space-time symmetric lattices with the temporal lattice size N t = 2, 



building upon a previous work fH)fl . Simulations have been carried out for a wide range 



of spatial lattice sizes, and finite-size scaling study of Lee- Yang zeros is used to find the 
location of the endpoint. We measure the Higgs and W boson masses around the endpoint 
and estimate the value of the Higgs boson mass at the endpoint. 

This paper is organized as follows. In Section |TJ we present the SU(2) gauge-Higgs model 
lattice action and outline our strategy for finding the endpoint through Lee- Yang zeros. In 
Section []TJ, following a brief discussion of susceptibility analysis, Lee- Yang zeros are exam- 
ined. Another approach to find the endpoint using the Binder cumulant is also described. 
In Section [IV| we present results of the zero-temperature mass measurement. Together with 
our result for the scalar self-coupling constant at the endpoint obtained through Lee- Yang 
zero analysis, this leads to the value of the Higgs boson mass at the endpoint. Sec. [V] is 
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devoted to conclusions. 



II. THEORY AND SIMULATION 

We work with the standard SU(2) gauge-Higgs model action given by 



s = J2 



E §th^ + E 2 ^ - pI - KpI - 1) 



(1) 



= iTr($t^$ x+A ), p 2 s Ixr^t^), ( 2 ) 

where U x ^ v is the product of link operators around a plaquette, (3 is related to the tree-level 
gauge coupling as (3 = 4/g 2 , k represents the Higgs field hopping parameter and A is the 
scalar self-coupling. We put the system on a space-time isotropic lattice of a size N t x iV 3 . 

Finding the endpoint of the first-order finite-temperature phase transition of the model 
requires finite-size scaling analyses to quantitatively distinguish the case of a first-order 
transition from that of a crossover as the coupling parameters of the model are varied. 



As the main tool, we employ finite-size scaling analysis of Lee- Yang zeros [[jl||T9[] on the 



complex k plane for fixed j3 and A @,|,|I7]. For a first-order phase transition, the infinite 
volume limit of the zeros pinches the real n axis, while they stay away from it if there is 
no phase transition. We also supplement this method with analyses of susceptibility and 
Binder cumulant. 

Our finite-temperature simulations are carried out for the temporal lattice size N t = 2. 
For the spatial lattice size we take iV 3 = 20 3 , 24 3 , 32 3 , 40 3 , 50 3 and 60 3 . The gauge 
coupling is fixed at (3 = 8. For the scalar self-coupling we choose five values, A = 
0.00075,0.001,0.00135,0.00145 and 0.0017235, which covers the range of zero-temperature 
Higgs boson mass 57 ^ Mh ~ 85GeV [fL5]|. For each value of A the scalar hopping parameter 



k is tuned to the vicinity of the pseudo critical point estimated by the peak position of the 
susceptibility of the Higgs field length squared p 2 . 

The updating algorithm is a combination of over-relaxation and heatbath methods jI3 



with the ratio of the two for the scalar part and the gauge part as specified in Ref. |L5| . 
We make at least 10 5 iterations of this hybrid over-relaxation algorithm at each coupling 
parameter point for each lattice size. The list of coupling values and statistics we use in our 
finite-temperature simulations are listed in Table |. 

We also carry out zero-temperature simulations to measure the masses of Higgs and W 
bosons around the endpoint of the first-order phase transition. For these runs an improved 
algorithm of Ref. |2(| is employed. Details of the runs and results are discussed in Sec. [TV} 

III. FINITE-TEMPERATURE RESULTS 
A. Susceptibility 

Let us first look at the susceptibility of squared Higgs length, 
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= V (<p 2 > - (p) 2 ) , (3) 

where V = Nf. The maximum value of the susceptibility at its peak, calculated by the 
standard reweighting technique []21| as a function of k, is plotted in Figure [I] against the 
spatial volume normalized by the critical temperature VT{? = N^/Nf. Errors are estimated 
by the jackknife procedure with the bin size of 10 3 -10 4 sweeps, which is listed in Table |. 

The slope for the smallest scalar coupling A = 0.00075 approaches unity for large volumes, 
which is consistent with a first-order transition, while that for the largest coupling A = 
0.0017235 tends to a constant, showing an absence of a phase transition. A continuous 
decrease of the slope for the intermediate values of A indicates that the endpoint of the first- 
order transition is located in between the two extreme values. Our range of spatial volumes, 
unfortunately, is not sufficient to pin down the critical value of A from the susceptibility 
data. 



B. Lee- Yang Zeros 

The determination of the endpoint of the finite temperature phase transition of the 
model, thus a characteristic feature of the phase diagram, is made by the use of the Lee- 
Yang zeros of the partition function Z fl8|IT9[| . Near the first-order phase transition point 
the partition function reads 

Z = Z s + Z b ex exp(-Vf s ) + eM-Vfb) , (4) 

where the indices s(b) refer to the symmetric (Higgs) phase and / stands for the free-energy 
densities. Near the phase transition point we also have 

fb = fs + Qf(« - K c ) , (5) 

since the free-energy density is continuous. One then obtains 

Z oc ex P [-V(f s + f b )/2) cosh[-Va(K - « c )/2] (6) 

which shows that for complex k Z vanishes at 

Im(re) = 2tt • (n - 1/2) /{Vol) (7) 

for integer n. In case a first-order phase transition is present, these Lee- Yang zeros move 
to the real axis as the volume goes to infinity. If a phase transition is absent the Lee- Yang 
zeros stay away from the real k axis. Thus the way the Lee- Yang zeros move in this limit is 
a good indicator for the presence or absence of a first-order phase transition. 

Calculation of the partition function for complex values of k is made with the reweighting 



method |2T] in both imaginary and real directions of k. In those cases where we have two 
ensembles with the same value of A and N s , but different k, we combine the two runs by 
setting the magnitude of the two partition functions to be equal at the midpoint between 
the two k's. 

In Fig. [| we show the absolute value of the partition function normalized by its value at 
the real axis on the complex k plane, 
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Znorm(^) — 



Z(KeK, Iran) 



Z(Rbk, 0) 



(8) 



for A = 0.00075 and N s = 60. The contour line of this figure is shown in Figure [| We 
observe three zeros in this case, whose distance from the real axis is roughly in the ratio 
1 : 3 : 5 as expected from ([/]) for a first-order transition. 

Let us call the zero nearest to the real axis as first zero, and denote its location by Kq. 
We search for the first zero by the Newton-Raphson method applied to the equation 

Z(Ren, Imre) = 0, (9) 

starting with an initial guess for kq obtained from the contour plot of Z norm (n). The error 
of kq is estimated by the jackknife method with a bin size given in Table |, i.e., the zero 
search is repeated for the set of partition functions calculated from each jackknife sample of 
configurations, and the jackknife formula is applied to the set of k,q. The results for n are 
given in Table |. We show in Fig. |] values of the imaginary part of the first zero ImKo(V) 
as a function of inverse volume. 

Finite-size scaling theory predicts that the volume dependence of the imaginary part of 
the first zero is given by a scaling form, 

ImKo(V) = «§ + CV~ V . (10) 

For a first-order phase transition, the infinite volume limit vanishes, Kq = 0, and the exponent 
takes the value v = 1. In the absence of a phase transition, k,q ^ and the value of the 
exponent is generally unknown. 

In Fig. |5| we plot results for clS cL function of A obtained by fitting the volume depen- 
dence of the first zero by the form ([!(]) (see Fig. [| for fit lines). Both Kq and v are taken 
as fit parameters, and the entire set of volume = 20 3 — 60 3 is employed. Filled symbols 
mean that they are directly obtained from the simulations carried out at the corresponding 
values of A. The points plotted with open symbols are obtained from the first zero of the 
partition function calculated by reweighting the partition function measured at the point 
where Kq with the filled symbol of the same shape is shown. The agreement of open sym- 
bols of different shapes within errors shows that reweighting from different values of A gives 
consistent results between the measured points. 

At small couplings A < 0.001, k,q is consistent with zero, which agrees with the result of 
Ref. [15| that the transition is of first order in this region. At large couplings A > 0.0013, 



Kq no longer vanishes, and hence there is no phase transition. In order to determine the 
endpoint of the phase transition, we take the three filled points at A = 0.00135, 0.00145 and 
0.0017235 directly obtained from independent simulations without A-reweighting, and make 
a fit with a function linear in A. This gives the position of the endpoint to be 

A c = 0.00116(16). (11) 



In Figure |6| we show the exponent of scaling function (10). The meanings of symbols are 
the same as in Figure [S]. For A > A c , where there is no phase transition, the exponent takes 
a value v ~ 0.75. Below the endpoint A < A c , the exponent shows some trend of increase, 
but not quite to the value v — 1 expected for a first-order transition. We think that this is 
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due to insufficient volume sizes used in our simulation, for which corrections to the leading 
1/V behavior are not negligible. 

To check this point we make an alternative fit of results for the first zero adopting a 
quadratic ansatz in volume given by 

lmK (V) = k c + CV- 1 + DV~ 2 } (12) 

and show the results for Kq in Figure [F| Clearly the infinite volume limit Kq starts to deviate 
from zero around A ~ 0.001, which is consistent with the estimate of A c above, albeit located 
at the lower end of the one standard deviation error band. 

We note that the quadratic ansatz ([12]) , formally the first three terms of a Laurent series, 
is expected to be correct in case of a first-order phase transition, for which ([!]) describes the 
thermodynamic limit. However, it is not a valid assumption in the region of A where there is 
no phase transition. Therefore, unlike the case of Fig. extrapolating the results of Fig. [7] 
from large to small values of A to estimate the location of the endpoint A c is not justified. 



C. Binder Cumulant 

Let us consider the Binder cumulant (cf. |p2]| ) of the space-like link operator, 

(L 4 ) 1 

BM = 1 " ^yip L s = _ J^L, >M (13) 

The infinite volume limit of the minimum of this quantity should deviate from 2/3 for a 
first-order phase transition, while it should converge to 2/3 beyond the endpoint. 

We evaluate the minimum of the cumulant as a function of k for a given A and volume 
using reweighting. We then use a scaling ansatz, 

= Bl + CV-\ (14) 

to extract the infinite- volume value B\ . 

In Fig. § we show — {B c Ls — 2/3) as a function of A, where the meanings of symbols are 
the same as in Fig. A change of behavior from non-vanishing values to those consistent 
with zero at A ~ 0.001 shows that the first-order phase transition terminates around this 
value. Linearly extrapolating the two independent data at A = 0.00075 and 0.001 yields 
A c = 0.00102(3) for the endpoint, which is consistent with the result ([11]) from our study of 
Lee- Yang zeros. Note, however, that only two measured points are available for the linear 
extrapolation. Therefore we can not make a statement on the goodness of the fit. For this 
reason, we conservatively take the Lee- Yang value (|IID as our best estimate of the endpoint. 



IV. CRITICAL HIGGS BOSON MASS 

To determine the physical parameters characterizing the endpoint, namely the ratio of 
the Higgs boson mass to the W boson mass and the renormalized gauge coupling g^, we have 



to perform zero-temperature simulations. As in Refs. |12|-|l4"|, we extract the Higgs boson 
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mass mu in lattice units from correlators of p x and L x ^. The W boson mass in lattice units 
mw is obtained from the correlator of the composite link fields 



W x = lMT r a+ + -U xk a x ), (15) 

r,k=l 

where r r is the Pauli matrix and a x is the angle part of $ x such that $> x = p x a x with 
a x E SU(2). 

Masses are extracted from the correlators fitting to a hyperbolic cosine plus a constant 
function. Simple uncorrelated least-square fits and correlated fits with eigenvalue smoothing 



proposed by Michael and McKerrell [[23[] are used. The application of this method is discussed 
in detail in Ref. [Q . 

The actual procedure of extracting the mass parameters is the following. First we deter- 
mine the reasonable time intervals for fitting the correlator data. The guideline is to choose 
as large an interval as possible with reasonable x 2 /d.o.f. value. For this purpose correlated 
fits with eigenvalue smearing are used. We find this to be necessary since the data are 
strongly correlated for different time distances. Having fixed the fitting time interval, we 
next carry out uncorrelated fits. To perform this fit, we divide the data sample into sub- 
samples, and estimate the errors of correlators from the statistical fluctuations of subsample 
averages. 

The best fit value of the masses is taken to be the number given by the uncorrelated 
fit. The value of the Higgs boson mass is obtained by fitting to a linear combination of 
the two different correlators for p x and L X}il . The errors on the masses are determined 
by jackknife analyses over subsamples. The masses obtained by the correlated fits with 
eigenvalue smearing are in all cases well within the error bars of the uncorrelated fits. 

Our zero-temperature simulations are carried out at two points given by (A, k = 
K c (X,N t = 2)) for A = 0.0011 and 0.00125 employing several lattice sizes to examine finite- 
volume effects. The run parameters and results for masses are collected in Table [TJ. The 
size of subsamples is typically 500 sweeps. 

Our results do not show significant volume dependence (see Fig. except for the two 
smallest spatial volumes N s = 8 3 , 10 3 for which somewhat different values are obtained 
compared to those of other volumes. We then discard those results and take an average over 



the rest of the volumes. This yields the values given in Table |TT[ Setting Mw = 80 GeV, 
we obtain 

M H = 70.9 ± l.lGeV (A = 0.0011) (16) 

M H = 76.8 ± l.lGeV (A = 0.00125). (17) 

Making a linear interpolation to the critical value A c = 0.00116(16) from the Lee- Yang zero 
analysis, we find 

M H<C = 73.3 ± 6.4GeV, (18) 

where the error is dominated by that of A c . 

From measurements of Wilson loops we also determine the values of the renormalized 



gauge coupling gji using the method described in Refs. |12]--|14f . The potential as a function 
of the distance R is fitted by 
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V(R) = -^e- MR + C + DG(M,R,L s ), (19) 
B 



where G(M, R, L s ) stands for lattice artifacts (cf. fll3fl). The potential is determined from the 
rectangular Wilson loops by fitting the time dependence with three exponentials. A stable 
fit is obtained in all cases. The potential is then fitted by ([19]) using all R values. Our results 
for the fit parameters and g\ for various spatial size lattices are shown in Table 0. We see 
that gn is constant within errors. The averaged values are given in Table |T|. The values 
do agree within errors, showing that our simulations for the two A values correspond to 
the same renormalized gauge coupling. Therefore the linear extrapolation to A c mentioned 
above is justified, since we use Higgs masses at equal renormalized gauge couplings. 

Finally, let us try to estimate the effect of fermions and the U(l) gauge boson on our 
result. We make this estimation through the perturbative expression for the parameter 
x = ^3/93 °f the dimensionally reduced model in terms of the physical parameters of the 
Standard Model p4|. Using our results for the Hiees boson mass and the renormalized gauge 



coupling, we find x c = 0.121 ± 0.020 for the endpoint. Including the effect of fermions and 
the U(l) gauge boson, this value corresponds to Mh, c = 80 ± 7 GeV. 



V. CONCLUSIONS 

We have studied the endpoint of the finite-temperature first-order transition of the SU(2) 
gauge-Higgs model on a space-time isotropic lattice of a temporal extension N t = 2. The 
results from Lee- Yang zero and Binder cumulant analyses show that the first-order phase 
transition terminates at A c = 0.00116(16) and turns into a smooth crossover for A > A c . 

Setting Mw = 80 GeV our result for the critical Higgs boson mass is M# c = 73.3 ± 6.4 



GeV. This is consistent within error with the value Mh, c — 74.6 ± 0.9 GeV [17] obtained in 
a 4-dimensional anisotropic lattice simulation for the same temporal size. The same work 
also reported that the critical mass decreases for larger temporal size, and extrapolates 
to Mh jC — 66.5 ± 1.4 GeV in the continuum limit. This value is consistent with the 3- 
dimensional result 66.2 GeV [[J. Thus results from various methods, in three and four 
dimensions, agree well. 



For a comparison with the experimental lower bound Mh > 87.9 GeV |55] for the 
Higgs boson mass, we need to include the effect of fermions and U(l) gauge boson. The 
good agreement of critical mass from the four- and three-dimensional simulations noted 
above imply that this may be made perturbatively, with which we find Mh, c = 80 ± 7 
GeV for our N t = 2 simulation. This value is about 10% larger, albeit with a comparable 
error, than the result M# jC = 72.4 ± 1.7 GeV in the continuum limit obtained from a 4- 
dimensional anisotropic study |T7] , possibly due to scaling violations. We also note that the 
3-dimensional approach reported the values M H c = 72.4 ± 0.9 GeV and M# )C = 72 ± 2 



GeV |10| . Combining all the available results, we conclude that the electroweak baryogenesis 



within the Minimal Standard Model is excluded. 
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TABLES 



TABLE I. Run parameters of finite temperature simulation and results of first Lee- Yang zero. 
Data used for analysis of susceptibility and Binder cumulant are marked with x and B, respectively, 
in the last column. 









(xlO 3 


sweep) 








A 


N s 


re 


iteration 


bin size 


Rereo 


Imreo 


use 


0.00075 


20 


0.129114 


100 


2 


0.1291133(23) 


0.0000477(20) 


X,B 




24 


0.129103 


100 


2 


0.1291068(12) 


0.0000285(11) 


X,B 




32 


0.129102 


100 


4 


0.12910273(91) 


0.00001351(45) 


X,B 




40 


0.129100 


100 


6.25 


0.12910086(72) 


0.00000762(26) 


A 5 




50 


0.129100 


120 


10 


0.12910041(51) 


0.00000411(17) 


X,B 




60 


0.129100 


180 


10 


0.129100308(303) 


0.000002321(51) 


X,B 


0.001 


20 
20 


0.129340 
0.129350 


100 
100 


2 
2 


0.1293472(15) 


0.0000605(16) 


B 

X 




24 


0.129330 


100 


2 


0.1293357(18) 


0.0000432(18) 


X,B 




32 


0.129328 


100 


2 


0.12933093(122) 


0.00002136(75) 


X,B 




40 


0.129327 


100 


2.5 


0.12932802(80) 


0.00001223(44) 


X,B 




50 
50 


0.129327 
0.129328 


100 
100 


4 
4 


0.12932797(37) 


0.00000763(26) 


B 

X 




60 


0.1293275 


180 


7.5 


0.12932743(37) 


0.00000489(18) 


X,B 


0.00135 


20 


0.129660 


100 


1 


0.1296888(34) 


0.0001167(36) 


X,B 




24 


0.129650 


100 


1 


0.1296619(29) 


0.0000819(32) 


X,B 




32 


0.129644 


100 


1 


0.1296465(20) 


0.0000542(20) 


X,B 




40 


0.129640 


100 


2 


0.1296426(15) 


0.0000293(11) 


X,B 




50 


0.129639 


120 


2.5 


0.12963782(137) 


0.00002016(88) 


X,B 




60 


0.129637 


120 


4 


0.12963754(68) 


0.00001299(78) 


X,B 

A. 5 


0.00145 


20 


0.129748 


100 


1 


0.1297482(35) 


0.0000885(38) 


X,B 




24 


0.129736 


100 


1 


0.1297384(20) 


0.0000567(20) 


X,B 




32 


0.129728 


100 


2 


0.1297318(15) 


0.0000328(12) 


X,B 






D 199794 


1 nn 

1UU 


9 

Li 


D 1 99797^1 (\ 1 51 

\J. ±Li<J I jLi i Ull lllj 1 


n nnnn9i 7i (qq) 

U.UUUU^ 11 11 ft? 1 


A) D 




50 


0.129722 


120 


2 


0.12972654(80) 


0.00001529(79) 


X,B 




60 


0.129724 


120 


4 


0.12972517(61) 


0.00001146(79) 


X,B 


0.0017235 


20 
20 


0.129980 
0.129990 


100 
100 


1 
1 


0.1299875(20) 


0.0000951(19) 


X,B 




24 


0.129980 


100 


2.5 


0.1299755(24) 


0.0000604(21) 


X,B 




32 


0.129966 


100 


1 


0.1299654(15) 


0.0000383(12) 


X,B 




40 


0.129968 


100 


1 


0.1299663(15) 


0.0000276(14) 


X,B 




50 
50 


0.129965 
0.129966 


100 
100 


2 
2 


0.1299616(14) 


0.0000207(16) 


X,B 
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60 0.129962 120 4 0.12996122(71) 0.00001585(74) X ,B 



TABLE II. Run parameters of zero-temperature simulations and results for masses in lattice 
units. 



A 




(xlO 3 sweep) 






K 


N'3 x N t 


iteration 


m H 


my/ 




8 3 x 20 


60 


0.2938(44) 


0.3583(41) 




in3 v o/i 
1U X Z4 


I 


n r >P,R r >( r >A\ 


U.oooU^oo J 


0.0011 


12 3 x 28 


49 


0.2844(46) 


0.3171(68) 


0.129416 


14 3 x 32 


34 


0.2838(34) 


0.3191(69) 




16 3 x 36 


26 


0.2851(62) 


0.3152(133) 




18 3 x 36 


26 


0.2887(47) 


0.3321(100) 




8 3 x 20 


60 


0.2806(42) 


0.3285(95) 




10 3 x 24 


75 


0.2764(33) 


0.3291(30) 


0.00125 


12 3 x 28 


49 


0.2884(38) 


0.2y92(51) 


0.129532 


14 3 x 32 


34 


0.2851(56) 


3037f58) 




16 3 x 36 


29 


0.2863(91) 


0.2941(64) 




18 3 x 36 


31.5 


0.2892(54) 


0.2965(71) 


TABLE III. Averaged masses in 


lattice units and renormalized gauge couplings 


from results in 


Table [D] excluding those for the two smallest volumes. 






A 


m H 


my/ 


Rhw 




0.0011 


0.2852(22) 


0.3202(41) 


0.8864(136) 


0.5712(27) 


0.00125 


0.2877(26) 


0.2988(30) 


0.9607(134) 


0.5768(33) 


TABLE IV. Summary of the fit parameters for the static potential and the renormalized gauge 


coupling. 










A 

K 


N 3 x Nt A 


M 


D C 


g\ = f ,A 




12 3 x 28 0.03495(58) 


0.3021(62) 


0.03941(68) 0.0968(2) 


0.5856(97) 


0.0011 


14 3 x 32 0.03435(52) 


0.2783(90) 


0.03673(45) 0.09672(21) 


0.5755(87) 


0.129416 


16 3 x 36 0.03406(30) 


0.2898(107) 


0.03975(28) 0.09632(13) 


0.5707(50) 




18 3 x 36 0.03394(22) 


0.2791(42) 


0.04061(262) 0.09633(3) 


0.5687(37) 




12 3 x 28 0.03561(46) 


0.2788(121) 


0.02814(310) 0.09751(31) 


0.5966(77) 


0.00125 


14 3 x 32 0.03456(70) 


0.2573(113) 


0.0353(57) 0.09766(39) 


0.5791(117) 


0.129532 


16 3 x 36 0.03386(31) 


0.2559(79) 


0.0416(35) 0.09740(12) 


0.5673(52) 




18 3 x 36 0.034442(35) 


0.2676(36) 


0.03831(41) 0.09704(4) 


0.5770(59) 
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FIG. 1. Peak height of susceptibility of p 2 against inverse volume normalized by critical tem- 
perature VTj} = Ng/Nf. Dotted lines are guides for eyes. 




FIG. 2. Absolute value of normalized partition function as a function of complex k for 
A = 0.00075 and N s = 60. 
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FIG. 3. Contour plot of Fig. |. 
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FIG. 4. Imaginary part of first Lee- Yang zero as a function of inverse volume normalized by 
the critical temperature. Solid lines are least x 2 fits with lm.Ko(V) = Kq + CV~ V . 
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FIG. 5. Imaginary part of first Lee- Yang zero at infinite- volume limit as a function of 
Higgs self coupling. Filled symbols are calculated without A-reweighting, while open sym- 
bols with A-reweighting from the filled symbol with same shape. Solid line is a linear fit to 
A = 0.00135,0.00145 and 0.0017235 (filled symbols). 
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FIG. 6. Exponent v of finite size scaling of first Lee- Yang zero. 
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FIG. 7. Same as in Fig. ||. Quadratic polynomial is used for fit instead of power function. 
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FIG. 8. Minimum value of Binder cumulant of L, at infinite volume limit as a function of A. 
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FIG. 9. Higgs and W masses in lattice units as a function of N s for A = 0.0011. 
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